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Abstract 

We present a semiclassical expansion of the smooth part of the 
density of states in potentials with some form of symmetry. The den- 
sity of states of each irreducible representation is separately evaluated 
using the Wigner transforms of the projection operators. For discrete 
symmetries the expansion yields a formally exact but asymptotic se- 
ries in h, while for the rotational SO(n) symmetries the expansion 
requires averaging over angular momentum as well as energy. A nu- 
merical example is given in two dimensions, in which we calculate the 
leading terms of the Weyl expansion as well as the leading periodic 
orbit contributions to the symmetry reduced level density. 
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1 Introduction 



In the analysis of quantum spectra it is common to decompose the full density 
of states into a smooth average term with a superimposed fluctuation, p(E) = 
p{E) + p osc (E). Each of these contributions has a distinct interpretation in 
terms of classical mechanics. The average or Weyl term is, to first order, 
understood as arising from orbits of zero length and is proportional to the 
allowed volume of classical phase space, while the oscillating term reflects 
the periodic orbit structure of the system. 

The smooth part of the density of states can be obtained semiclassically, 
by expanding the density operator p = S(E — H) as a series in h, 

p(E)=trp = J2A n (E)h n (1) 

n 

which we shall refer to as the Weyl expansion. Formally this gives an exact 
expansion of p, and we should recover from this series both the average and 
the oscillating parts of the level density. On the other hand, any truncation of 
Eq. ([I]) is only a slowly varying function of the energy and does not contribute 
an oscillating term to the level density. The resolution is that in general we 
expect the Weyl expansion to be asymptotic in % and it has been argued [|l| 
that information about the oscillating part of the density is encoded in the 
divergent part of the series. Indeed, by resumming the remainder of the series 
in a billiard system oscillatory contributions given by the shortest periodic 
orbits are found. 

The purpose of this paper is to expand the Weyl term in the presence 
of either a discrete or a continuous symmetry. Such a symmetry implies 
conserved quantum numbers so that the density decomposes into a sum of 
partial densities, p(E) = J2Pa(E), where a labels a complete or partial list 
of conserved quantum numbers. We are more often than not interested in 
the reduced densities p a (E), rather than the full density of states. For exam- 
ple, a quantum system composed of identical particles (e.g. fermions) have 
wavefunctions that are totally odd or even under particle exchange and the 
density of states is given by the density of the antisymmetric or symmet- 
ric representations of the symmetric group. Another example would be an 
analysis of an experiment in the presence of selection rules. For instance, 
we might have an experimental probe of a spherically symmetric potential 
which is only sensitive to states with a specific value of J, so we would need 
to know the density of such states. 
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A practical use of the Weyl expansion is in measuring the reliability of 
both experiments and numerical calculations. By comparing the reported 
density of states to the Weyl expansion one can decide if one has, for instance, 
missed levels. In the presence of symmetry, it is easier to do experiments and 
numerical calculations one symmetry representation at a time and we need 
the symmetry reduced Weyl expansion for comparison. Finally, it is common 
to look at various statistical aspects of spectra such as spacing distributions 
[0] and to do so requires that we first unfold the spectrum by removing 
the Weyl term from the density of states. Since these statistics are only 
meaningful when computed separately for each symmetry representation, it 
is necessary to know the symmetry decomposed Weyl expansion. 

As an asymptotic series, the Weyl expansion is most accurate when trun- 
cated one term before its least term and therefore the number of terms which 
we are allowed to include in the smooth part depends on %. We will not go 
into any details on resumming the Weyl expansion. Rather we adopt the 
practical point of view that we can simply add the leading oscillatory part 
of the density to the smooth density of states. Thus, for the partial density 
p a of states belonging to the irreducible representation a of the symmetry 
group we write 

p a {E) « p a {E) + (E) (2) 

where the average term is given by the truncated Weyl expansion. The 
oscillating term p a:OSC has previously been treated for both discrete [3-5] and 
continuous || |7| symmetries. To first order it yields a trace formula similar to 
the Gutzwiller trace formula but with the periodic orbits given in a reduced 
phase-space. 

The Weyl expansion in the presence of discrete symmetries was previously 
discussed in Refs. [8-10]. Our formalism is the same as that of Refs. || [5j. In 
Ref. || however, the discussion was limited to a single or several commuting 
reflections while Ref. [0] was specific to the symmetric group for identical 



particles. In Ref. |T0[ the first few terms of the expansion was found for 
billiards. Here however, the discussion was specific to quantum billiards and 
it does not seem easy to generalise it to potential systems. 

The structure of this paper is as follows. In section 2 we discuss the 
mechanics of the decomposition of the level density into the irreducible rep- 
resentations. Both the density operator and the projection operators are 
evaluated using the Wigner representations of the operators so this topic is 
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briefly reviewed. We discuss discrete groups in section 3 and the continu- 
ous group SO(n) in section 4. All point group elements can be expressed 
as a multiple reflections, in n dimensions there are at most n independent 
reflections. The reduction of more complicated group elements to multiple 
reflections allows for an ^-expansion of the operators. For continuous groups, 
we find that the semiclassical limit and the limit of small rotation angles are 
noncommuting. The results we obtain are found to only be sensible if we 
further average over angular momentum. In section 5 we discuss the Henon- 
Heiles potential which has a C^ v discrete symmetry. We compare the results 
of sections 2 and 3 with the actual spectra. By studying specific combinations 
of the reduced level densities, we can isolate the contribution from different 
group elements. We find that we can also discern the effects of the shortest 
periodic orbits on the oscillating density of states, as discussed in section 6. 
Section 7 contains some concluding remarks. 

2 The Formalism 

We are interested in the general class of quantum Hamiltonians, 

for which the potential V(q) is invariant with respect to the symmetry trans- 
formations of a group G. In the present paper, we consider all groups G that 
are subgroups of the group of orthogonal transformations in n dimensions, 
0(n). The groups G can include improper rotations, and can be either dis- 
crete or continuous. These groups clearly cover many physical applications, 
they include the crystallographic point groups and the symmetric group S n 
of n identical particles. 

The Hamiltonian commutes with all the group operators of G and the 
eigenstates of the Hamiltonian can be classified according to the irreducible 
representations (irreps) of the group. This means that we can decompose the 
full density of states into densities of states belonging to each irrep, 

p(e) = J2p«(e), (4) 

a 

where the index a runs over all the irreps of G. Corresponding formulae 
apply for p(E) and p osc (E). Each partial density p a (E) may be obtained by 
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projection, 
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= J dfi(g)xa(9)Ul 1 1 dfx(g). (6) 



The first expression applies when the group G is discrete; d Q is the dimension 
of the irrep and \G\ is the order of the group. The sum is over all of the group 
elements g G G and x«(<?) is the character of element g in representation a. 
The operator U g transforms wavefunctions as prescribed by the group element 
g. The second expression applies when G is continuous, in which case the 
sum is replaced by the integral over the group elements with the invariant 
Haar measure dfi(g). 

In the semi classical limit, average densities are conveniently described by 
phase space integrals through the use of the Wigner representations of the 
various operators. The Wigner transform of a generic quantum operator A 
is defined as 



This has the appealing property of being a complete representation of the 
quantum operator defined in the classical phase space. The phase space is 
assumed to be 2n-dimensional so q, p and x are all n-dimensional vectors. 
Traces have a simple evaluation in this representation since 



These equations do not generalise to more than two noncommuting operators 
as can be seen from the fact that the ordering of the operators in the traces 
is important while the ordering in the integrals is irrelevant. Instead we use 
the product rule, 




(7) 






(9) 
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where the Moyal operator A is given by 

A=V g ■ V P - V P • V q (10) 

and the arrows indicate which way the derivatives should be applied. Eqs. (|8|), 
(HD and (|H]) allow us in principle to evaluate the trace of any operator of the 
form ABC ■ ■ ■ given the Wigner transforms of the individual operators. A 
final property of the Wigner transform which we will use is that Hermitian 
conjugation transforms into complex conjugation, (A^)vk = (Aw)*- In par- 
ticular, this implies that the Wigner transform of any Hermitian operator is 
real. 

Using the definition (0) we find that the Wigner transform of the quantum 
Hamiltonian (|3|) is simply the classical Hamiltonian 

*Mq,p) = y + ^(q)- (ii) 

To evaluate the partial densities of Eq. ([5]) we need the Wigner transform of 
the density operator 8(E — H)w- However, it is not true that 5(E — H)w 
equals 5(E — Hw). The difference between these expressions gives the higher 
order contributions to the full density of states. These have been worked out 
in a particularly elegant manner in Ref. The density operator can be 
written as the inverse Laplace transform, 

where the Wigner transform of the evolution operator is expanded in h, 



e 



-j3Hw 



1 /-fr 2Xm 



Sm t AMl (13) 



The coefficients A m are simple polynomials in j3 and obey the recursion 
relation, 



d P \ 2k J [ 



A = 1, (14) 
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from which one may easily generate the series. With the Hamiltonian (3), 
the first few terms of the density operator are 



5{E-H) W = 5{E-H w )-%-V 2 V^-5{E-fr 



w) 



dE 2 

+- ((p ■ V)V + (W) 2 ) -^6{E - H w ) + 0(ft 4 )(15) 
a result already obtained by Wigner in 1932 fl2]| . An alternative approach, 



which yields the same expansion, was developed in Ref. 13|| . In calculating 
the density of states, the derivatives with respect to energy can be taken 
outside the phase space integral. 

The semiclassical expansion of the density operator only contains even 
powers of h. This is in contrast to the expansion for quantum billiards, 
which in general contains all integer powers of h, the first odd term arising 
from the Dirichlet or Neumann boundary conditions. In general, quantum 
billiards require a separate treatment; the semiclassical expansions of the 
Wigner transforms above do not converge when the billiard boundaries are 
approximated by steep but smooth potential walls. We will find that the 
projection operators @ also contain odd powers of h in their Wigner expan- 
sions. In simple cases, the various irreps may be realized by internal hard 
wall boundaries thus explaining the odd powers. However, the semiclassical 
expansion as stated does not allow for singular boundaries. Instead, we will 
directly evaluate the Wigner transforms of the projection operators P a or 
equivalently of the transformation operators C/j. 

Eq. (p~5|) appears to contain only the smooth part of the density of states. 
The derivation is certainly exact so we should expect the trace of the density 
operator to contain the oscillating term in addition to the average term. 
The resolution of this apparent paradox lies in the concept of resurgence [p]] . 
The Weyl series is an asymptotic expansion in h and information about the 
oscillating density of states is contained in the very high order terms. In that 
sense, there is no loss of information in Eq. flT5|), only when we truncate the 
series we disregard the oscillations in the density of states. The truncated 
Weyl series gives the average density of states. It follows, that the separation 
between the average and oscillating terms of the density of states is somewhat 
arbitrary. In principle, we could keep all terms in the expansion until the 
terms start to grow. In practise, we are satisfied with a few terms, e.g. all 
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terms up to order 0(h n+1 ) in Eq. ( |T5D or equivalently, all terms up to order 
0(h) in the density of states. 

3 Discrete Groups 

The point groups in n dimensions are subgroups of the group of orthogonal 
transformations 0(n) and the group elements can be represented by n x n 
orthogonal matrices acting on the coordinates. The eigenvalues of an orthog- 
onal matrix R € 0(n) have unit norm and either come in complex conjugate 
pairs, A = exp(±£0) or are real, A = ±1. An eigenvalue A = — 1 corresponds 
to a reflection in a hyperplane of dimension n — 1 while the complex eigen- 
values exp(±z0) correspond to a rotation by the real angle in a hyperplane 
of dimension 2. It follows that by an orthogonal transformation, R can be 
brought into a block-diagonal form, 



(16) 



where the blocks are either one- dimensional = ±1 or two-dimensional of 
the form 

/COS0, -sm0A (17) 
V sin (pi cos (pi J 

Thus any orthogonal transformation in n dimensions can be constructed out 
of k successive, independent transformations, 

R = Rk ■ ■ ■ R2R1 (18) 

where each transformation Ri is either a reflection or a rotation. The trans- 
formations Ri commute since they act in orthogonal hyperplanes and the 







. 
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decomposition (|T8|) is unique up to an irrelevant ordering of the k transfor- 
mations. In three dimensions we infer that there are only four types of group 
elements, the identity J, a simple reflection R a , a proper rotation R^ and 
an improper rotation R a R<f,- For the improper rotation, R a is a reflection 
through the rotational plane. 
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Corresponding to each coordinate transformation Ri, there is a quantum 
operator Ui = U(Ri) which operates on the Hilbert space. These operators 
have the same multiplication properties as (|T8D and are defined by 



U^{v)=^{R- l v). (19) 

We can now evaluate the density of states belonging to the irrep a of a 
discrete symmetry group G. From Eqs. @, (|6|) and (^) we have 

P^ E ) = |E Xa(9) J 0l)wS(E - H) w . (20) 

The semiclassical expansion of 5(E — H)w is given by Eq. (|H>p. The Wigner 
transform (U g )w is most easily evaluated in a coordinate basis where the 
matrix representation of the group element g attains the block-diagonal form 
(P~6[) . By direct application of the definition, Eq. ([F]), it is clear that with the 
decomposition (|18D the Wigner transform {U g )w also factorises as 

(U g ) w = (U k )w-(U 2 ) w (U 1 ) w . (21) 

Each (Ui)\y{<i, p) is a function of only those coordinates and momenta that 
are not left invariant by the transformation. The factorization ( pT|) could also 
have been obtained by repeated applications of the product rule (|9|), using 
the orthogonality of the invariant hyperplanes of the transformations R{. 

It remains to evaluate the Wigner transform of the basic operators tj}. 
The identity element I is trivial since 

(Ui)w = 1, (22) 

corresponding to an empty product (|2~T|). From Eq. (|7|) it also follows that 
the Wigner transform of a reflection through the qi = hyperplane is given 
by 

(U a )w = /^ 1 e- i ^(g 1 + ^|-g 1 + ^) (23) 



7Th5(q 1 )8(p 



The quantum operator corresponding to a rotation by an angle <fi in the q\ — q<i 
plane is given by = exp (—i(f>L 3 /h). Its adjoint has the Wigner transform, 



e 2itan(</>/2)(<jip2-92Pl)A 
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The Wigner transforms of the identity and reflection operators automatically 
lend themselves to a semiclassical expansion. The rotation operator on the 
other hand, poses a problem as its Wigner transform (]23j) contains an essen- 
tial singularity in h = 0. This may however, be circumvented in the following 
manner. 

Each two-dimensional rotation by an angle can be constructed from 
two consecutive reflections about two axes (located in the two-dimensional 
hyperplane) that intersect at an angle 0/2. Two such reflections do not 
commute, the opposite ordering would amount to an inverse rotation. From 
the product rule (|9]) and the Wigner transform of a reflection operator (B3) 
we may write 

{U\) w = (7rh) 2 8(q' 1 )5{p[)e- in ^ 2 5(q 1 )6(p 1 ) (25) 

where q\ = and q[ = denote the non-orthogonal invariant hyperplanes 
of the two reflections. The essential singularity is now only implicit through 
the derivatives of the distribution functions 5(q), 5(p), and we shall simply 
expand the expression ( ^5|) in H. Since the q[ = hyperplane is rotated from 
the qi = hyperplane by the angle 0/2 we obtain, 

{U\) w = (7r^) 2 5(cos(^) gi + sin(^)g 2 )5(cos(^) Pl + sm(^)p 2 )e-^ 2 5( qi )5(Pi) 
U,1)2 -Mq2)S( P2 ) 




d d d d 



— m: \ £. * j ydq 2 dp 1 dp 2 dq l 
To leading order in H we have derived the simple result 



5( qi )S( Pl ). (26) 



« .\! >/0 M q l )8(p 1 )5(q 2 )8(p 2 ) (27) 
v sm (0/2) 



also valid for billiard systems [jT0| . The higher order terms containing the 
partial derivatives may evaluated when we perform the phase-space integral 
of the reduced density of states, Eq. (|20|). That is, integrating by parts we 
obtain for the "rotational density" , 

dpdq 

" iw 



= tr{Ul5(E -H)} = J j^~ n 0l)w6(E - H) 
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/ (/p ' /q S{qi)8(pi)6{q 2 )5(p 2 ) 



sin 2 (0/2)7 (2ixhy 

00 1 fih \ m f d 2 d 2 \ m 

xE^U^ W 2)) g^-g^] HE-HM28) 



m=0 



All the derivatives now act on the Wigner transform of the density operator, 
8(E — H)w With the expansion dl5f) of the density operator, it is straight- 
forward, if tedious, to write down the integrals of Eq. (28) in closed form. 



We do not do so here but only discuss the leading order contributions to each 
partial density of states p a . 

The identity operator is present in any symmetry group G. In the semi- 
classical expansion, this gives the leading order term for all representations. 
Recalling that the character of the identity element is simply the dimension 
of the representation we get as a first approximation 

Pa(E) « ^p(E), (29) 



where p(E) is the total density of states. This was also noted in Ref. [p~0| 



where the author argued that this expression holds equally for potential 
systems and billiard problems. The next to leading order term is given by 
single reflection operators if they are present in the group. Rotation operators 
contribute to the second order corrections as seen in Eq. (|26|). As mentioned 
in the last section, reflection symmetries lead to the odd powers of H which 
are absent in the full density of states. 

In the general case, more complicated group elements in the sense of the 
factorisation (pip contribute to higher orders in h. That is, for each group 
element g, the lowest order in h of the expansion of (U g )w is given by the 
number of independent reflections contained in g, a rotation counting as 
two reflections. Therefore, the power of h with which U g contributes is the 
codimension of the set of points which is invariant with respect to the group 
element g. Loosely speaking, it is the invariant points which contribute to 
the zero length pseudo-orbits from which the Weyl term can be evaluated, 
as argued in Ref. [[HJ. For example, the points which are invariant under 
a finite rotation are those along the axis of rotation and comprise a set of 
codimension two. Indeed, the corresponding operator contributes with two 
added powers of % relative to the identity operator. Note in this case, that 
the limit of a small angle 0, the contribution diverges. This is because all 
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points become approximately invariant in this limit. The difficulty is that 
the limits h — > and — > do not commute; we discuss this problem in the 
next section. 

To conclude this section we give the explicit expansions for a two-dimensional 
potential system, keeping only terms with non-positive powers of H. The re- 
mainder of the smooth density of states then vanishes in the limit of h — > 0. 
Corresponding formula for higher dimensional systems may easily be ob- 
tained with the help of the factorization (|20|). It is convenient to work with 
the integrated density of states defined by 

N a (E) = [ E dE' Pa (E') (30) 

J — oo 

which is the number of states in representation a with energy less than E. 
We then find from Eq. (|20|) 



Na(m = ^iY,Xa(9)N g (E), (31) 

ICrl g 



where we have defined 

N g (E) = f jE'\ 0l)w5(E - H) w . (32) 

The quantity N g (E) is the contribution to the various integrated densities of 
states from group element g. We obtain up to order h, 



Nj(E) 



dpdq 

2 



(2vrft) 



e(E - h w ) - 1 -^- 2 e(E - %)W(q) 



h 2 d 3 



+ Ya~Ie* q{e ~ Hw) V (p ' v) V + (VF)2 



+ o(h 2 



N„{E) = \J^- &(E - H w )\ qi=pi=0 + 0(h) 
1 



N + {E) = 4sin 2 ^/2) Q(^-^)l q =p=o + Q(^) ; (33) 

where Q(x) is the Heaviside step function and the reflection a is defined to 
be in the q\ — axis. 
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4 Continuous Groups 



Much of the discussion of the discrete point groups also applies to the contin- 
uous point groups. We begin by studying the abelian group SO (2). However, 
as we argued in the last section, we may generalize the results to the group 
of rotations in n dimensions, SO(n), by use of decompositions of the form 
fll8p . The SO (2) group is also relevant to integrable Hamiltonians where the 
SO (2) symmetry is manifest in the action-angle variables since the Hamilto- 
nian does not depend on the angles. 

We define the rotation to be in the q± — g 2 plane and denote by L 3 the 
corresponding angular momentum operator. The irreps are then labelled by 
the angular momentum eigenvalues m and the projection operator is given 
by Eq. ©, 

P m = ^- I dd e- m9 Ul (34) 



2tt 

with the rotational operator 

U} = exp(i6L 3 /h). (35) 

From the form of the rotational operator it is clear that for continuous groups 
the semiclassical limit, % — > 0, does not commute with the small angle limit, 
9 — > 0. This is also obvious from the expansion of the "rotational density", 
Eq. (|27D, where all terms diverge in the limit of small angles. Instead, we 
need to project before we expand semiclassically. The Wigner transform of 
the rotation operator is given by Eq. (|24]). 

f a e 2itan (0/2)L :i /h 

U ^w= cos* (0/2) ' (36) 

where L 3 is the Wigner representation of L 3 and equals the classical angular 
momentum qiP2 — c[2Pi- The density of states of angular momentum (L 3 ) = m 
then reads 

-ki^h ^kr m{en)L " H {e -^- (3?) 

This expression may be brought into a more convenient form for semi- 
classical expansion by substituting the variable z = tan (8/2)/h to obtain 
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We separately consider m = and m^O. For m = the integral over z is 
trivial, so the density of zero angular momentum states becomes, 

Po(E) = J hd(L 3 )5(E - H) w . (39) 

This result has the simple interpretation that the m = states correspond 
directly to the L 3 = region of phase space. If we are only interested in 
such states we can ignore the rest of phase space and quantise directly on 
the L 3 = manifold. The result may be slightly surprising as it seems to 
neglect quantal fluctuations around that manifold. 

For non-zero values of m the integral can also be evaluated exactly but 
it will not lead to transparent results. Instead we expand the integrand in a 
formal series 

(l-iz%\ m ( ™(-izh) 2k \ 

exp —2imzh 



\l + izhj \ ~ %k + l ) 

= e~ 2mizh (l + 2 -^{-izhf + --^ (40) 

from which we obtain the desired semiclassical expansion of the density of 
states belonging to the irrep m, 



Pm{E) = J (%5{L 3 - mh) + ^5®(L 3 - mh) + ■ ■ •) 5(E - H) 



w ■ 



rf p rf q u s ,t , m 

(2irh) n V 

(41) 

As for m = we have an identification between states with a given value of 
m and the corresponding L 3 = mh region of phase space. However, there are 
corrections in the form of derivatives of delta functions. It is not enough to 
consider just the L 3 = mh manifold; we must also understand its neighbour- 
hood in order to evaluate the derivatives. In addition, it should be noted 
that the expansion in Eq. ( |4"0"| ) is convergent only if \ihz\ < 1 so strictly 
speaking it is incorrect to integrate z with this expansion all the way to in- 
finity. However, the error in doing so is small since the integral ( |3lf ) is highly 
oscillatory for large z. Furthermore, the oscillations are washed out if we 
further average over a range of m values. Therefore, we interpret Eq. (£y]) as 
the average density of states in both energy and angular momentum. Both 
of these averages are over ranges which are quantum mechanically large but 
classically small. Therefore this result only holds for m ^> 1. This may be 
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expected since it simply says that semiclassics works in the limit of large 
quantum numbers. 

The lowest order term of the expansion (IT) was obtained in Ref. H by 



integrating over the short orbits in a symmetry reduced Green function. The 
result here is slightly more general since we have shown, in principle, how to 
obtain the higher order terms. In addition, we have shown that the case of 
m = is special since the naive expectation is exact in this case. 

The generalisation to higher dimensional systems is now straightforward. 
If the Hamiltonian has the symmetry of k independent SO (2) groups, the 
eigenstates can be labelled by k quantum numbers. To leading order, the 
partial densities are 

Pm{E) = I (2^ " ~ (42) 

where m is the k dimensional vector of integers which labels the represen- 
tations of the k groups and L is the vector of the corresponding classical 
generators. 

Rotational invariance in three dimensions is described by the symmetry 
group £0(3). The irreps are labelled by the angular momentum I and are 
(21 + l)-fold degenerate. The characters of each irrep I is Xi{@) = sin((/ + 
\)6)/ sin(|) so that the projection operator becomes ]TT] 



/, = 2 2±1 [ da /^(g) sin «'+pW (43) 

An J J 2 sm(|) 

Here Q gives the direction of the rotation axis and Lq(9) is the corresponding 
angular momentum operator. The Wigner transform of the rotation operator 
is the same as in the 5*0(2) calculation and is given by Eq. 



2itan(0/2)L n /7i 

( iL n e/n\ _ f (aa) 

V )w cos 2 (^/2) ' { } 

where Lq now denote the value of classical angular momentum along the 
rotation axis. Substituting this for the rotation operator in Eq. (fi3|) and 
integrating over Q we obtain the Wigner transform of the projection operator, 

< fl »«' = — L ie L { i i am ( 2 tan( 2K > ) • < 45) 
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Here L is the length of the angular momentum vector. As in the two- 
dimensional case, we substitute for 9 the variable z = tan (9/2) /h and expand 
in h, 

2 (zh - ^(zhf + ■ ■ (46) 

To first order in h the Wigner transform of the projection operator takes the 
simple form, 

(P,V = <2f + i)^d + 

= (2( + l)fi 2 ,5(L' 2 - (( + i) 2 fi 2 ) (47) 

and the density of states of angular momentum I is given by 

Pl(E) = (21 + 1) / ^ h>5 (L 2 - (/ + 5(E - H) w . (48) 

For any L we have the result that we must evaluate quantities on the man- 
ifold L = (I + l/2)h. As before, the corrections to this expression involve 
derivatives evaluated on the manifold. Note that unlike in the case of SO(2), 
there is nothing special about the L = states. 

For higher dimensional groups SO(n), each orthogonal transformation 
decomposes as in Eq. (^) into a sequence of mutually orthogonal rotations 
Ri = exp(i9iL{) and we must determine the Wigner transform of a product 
of the form exp(^iLi) exp^^I^) • • •• This is a trivial extension of what 
we have done for 5*0(2) so in principle we have worked out the semiclassical 
expansion. Of course, for each group SO(n), one must include the characters 
specific to that group. 

5 A Numerical Example 

We shall verify the results of section 3 by evaluating the partial density of 
states in a simple symmetric potential system in two dimensions. Symme- 
try transformations in higher dimensional systems may, with respect to their 



6» = 2tan _1 (^) = 
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semiclassical expansions, be reduced to a sequence of reflections and rota- 
tions, already present in two dimensions. We use the Henon-Heiles potential 



V{x,y) = \{x 2 + y 2 ) + x 2 y-^y i 

= -r 2 + -r 3 sin39, (49) 
2 3 

first introduced as a model of celestial mechanics. The potential can be in- 
finitely negative but if E < 1/6, motion near the origin remains bounded and 
for small energies resembles that of a two dimensional harmonic oscillator. 

The potential is symmetric with respect to rotations by an angle = 
±2-7r/3 and with respect to reflections through three symmetry axes. To- 
gether with the identity element, these comprise the three classes of the 
six member symmetry group C% v . Note that this is isomorphic to the per- 
mutation group S3 which arises when we consider three identical particles 
interacting in one dimension. These groups have three irreducible represen- 
tations. Two are one dimensional and are labelled A\ and A 2 . The third is 
two dimensional and is labelled E. The character table is shown in Table 1. 

Using Eq. ([Jl]) we find 

N M {E) = l[N I (E)+2N (t> (E) + 3N (T (E)] 
o 

N M {E) = ^[N I (E) + 2N4E)-3N a (E)} 

N E (E) = 2 -[2N I {E)-2N^E)} (50) 
o 

We have used the fact that the contribution from each member of a class 
is identical. Note that the leading order contribution has the densities in 
the ratios 1:1:4 as we expect from Eq. (P9"D- We could calculate these three 
quantities separately in order to compare with the numerics. However, it is 
more instructive to isolate the contribution N g (E) from each symmetry class 
by inverting Eq. ( ^0|) as follows, 

Nj(E) = N Al (E) + N A2 (E) + N E (E) 
N a (E) = N Al (E) - N A2 (E) 

N^E) = N Al {E) + N A2 {E)-N E {E)/2. (51) 
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= rrf e [e-Wv+v -vn) (52) 



The first of these is just the total density of states and we do not consider it 
any further. The semiclassical expansions of the other two terms are given 
by equation (|3"3"D. For the Henon-Heiles potential they are to leading order 

O ft? f 1 ^ , 1 ..2 1 ..3- 

N(f>(E) = i (53) 

The reflection term N a (E) has to be be evaluated numerically. Note that 
the rotation function N^E) has no corrections at finite energy; all correction 
terms are of the form of derivatives of the stepfunction 0(E), which vanish 
identically for E > 0. 

We have calculated the exact quantum spectrum for % = 0.004 and eval- 
uated the smooth density of states by convolving the exact density of states 
of each representation with a Gaussian of width w so that 

p a (E) = - *rdE?e-l*-*ri** Pa {E?) 

V ll\W A J-oo 

where E n>a is the energy of the n'th state in representation a. We use the 
notation p for the smoothed density of states in order to distinguish it from 
the semiclassical average p. For consistency, we must smooth both p and 
the exact density of states; these smoothed quantities should agree when the 
smoothing width w is large. 

The results for the reasonably large value of w = .01 are shown if Fig. 1. 
The solid lines are evaluated from the exact spectrum while the dashed lines 
are the results shown in equations ( |52|) and fl5"3). In both cases the agreement 



is very good, in fact for N a (E) the dashed lines lies directly on the solid line 
and can not be discerned. Also note that the behaviour at small energy 
is modified by the smoothing; for N^(E) the curves begin at the value 1/6 
rather than 1/3 while for N a (E) the curves do not begin at the origin as they 
would for the unsmoothed case. There is also discrepancy at the high end 
of the spectrum due to the fact that the calculation of the quantum energies 
becomes unreliable near the threshold of E — 1/6. This affects N$ more 
than N a (E) since it relies on a higher order of cancellation (note the scale) 
and is therefore more sensitive to inaccuracies. 
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6 Effect of Periodic Orbits 



As shown in Fig. 1, the discussion of the previous sections accurately gives 
the smooth level density for each representation. To get this agreement it 
was necessary to convolve the delta functions of the exact spectrum with 
a Gaussian in order to smooth out the spectrum. However, as the width 
w of the Gaussian is decreased, we should expect first to see the periodic 
orbit structure in the spectrum and in the limit w — > to recover the exact 
spectrum. Here we consider only a relatively large smoothing width, so that 
only the shortest periodic orbits are apparent. 

We will discuss the periodic orbit structure not in the full configuration 
space but rather in one of its fundamental domains which tessellate the full 
configuration space. Dynamics in the fundamental domain are defined as the 
dynamics in the full domain mapped onto the fundamental domain. In par- 
ticular, periodic orbits of the full dynamics always map onto periodic orbits 
in the fundamental domain. Conversely, periodic orbits in the fundamental 
domain correspond to fractions of periodic orbits in the full domain. These 
fractional periodic orbits have endpoints which are related by one of the group 
elements. We can thus associate a group element to each periodic orbit in 
the fundamental domain. However, the group element depends on the choice 
of fundamental domain and is only specified up to its class. Orbits along the 
boundary of the fundamental domain are exceptional since group elements 
from more than one class may connect their endpoints in the full domain. 
The treatment of boundary orbits is intricate and, since we do not consider 
them in what follows, we forego an involved discussion. 

For interior orbits (i.e. non-boundary) the oscillating term of the partial 
densities is given by p QiOS c = — im {fi , a}/ 7r where [|J 

J oo oo 

9a{E) = -f E ^EE exp (ir [Sj/h + i(k + l/2)( Uj + ivj)]). 

ieppo r=i k=o 

(55) 

The index j labels all the primitive periodic orbits in the fundamental domain 
and the gj is the group element associated with orbit j. The index r counts 
the repetitions of the primitive orbits. T, 5*, u and v denote the period, 
action, stability exponent and angle respectively, and the bar denote the 
values when we consider the dynamics in the fundamental domain. This is a 
direct generalisation of the Gutzwiller-Voros formula [pi in the presence of 
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a discrete symmetry. 

Similarly to Eqs. floTp, we may define the "class greens functions", 



9i(E) 
9<p(E) 



g Al (E) + g M (E) + g E (E) 
g Al (E) - g M {E) 
g Al (E)+g A2 (E)-g E (E)/2. 



(56) 



which give the oscillatory contributions to the class densities. 

Like any bounded Hamiltonian system, the Henon-Heiles potential has an 
infinity of periodic orbits. However three stand out as special because their 
periods are relatively short. In Fig. 2 we show the configuration space with 
the three periodic orbits and the fundamental domain given by the bottom 
left triangle. In the figure, we also show the continuation of the periodic 
orbits in the full domain as dashed lines. Following Ref. [JHJ we label the 
three periodic orbits A, B and C. In the fundamental domain, orbit B begins 
at the short side, reflects off the hypotenuse to the long side where it turns 
around and retraces itself. Its period in the fundamental domain is half that 
of the full domain and we associate it with the class of reflections. Its double 
repetition has a full period and we associate it with the identity class. Orbit 
C starts at the short side, goes to the hypotenuse and is reflected back. It 
has one third the period of the full orbit and is associated with the rotation 
class. Its double repetition is also associated with this class but its triple 
repetition has the full period and is associated with the identity class. Orbit 
A runs down the short side to the centre and then up the hypotenuse. It then 
turns around and retraces itself so its period in the fundamental domain is 
the same as in the full domain. Since it is a boundary orbit it is associated 
with both the identity and the reflection classes. 

In the full domain, all three orbits have a period close to 2tt, the value 
they approach in the limit E —>■ 0. A and B are both triply degenerate while 
C is doubly degenerate because of its time reversed partner. However, in the 
fundamental domain all three orbits are non-degenerate. Orbit A fluctuates 
between being stable and unstable as the energy is varied. Orbit B is unstable 
for all energies and orbit C is stable for all energies up to E ~ 0.145 where 
it suffers a tangent bifurcation. In Fig. 3, we show the periods T as well as 
the stability angle v or exponent u for the two orbits B and C. The stability 
angle of the unstable orbit B is a constant v = Att, corresponding to a Maslov 
index 4. 



20 



The "rotational greens function" g<f>{E) probes the class of rotations and 
of the three short orbits only has contributions from C . Each third multiple 
modulus 1 or 2 of repetitions of the primitive orbit in the fundamental domain 
contributes to the density, and we obtain 

rp OO OO 

9*{E) = 4" £ £ {^ +r) + ^ (|+r) ) (57) 

where a = S/h — v (k + 1/2). Since long orbits will be damped by the 
smoothing parameter w, we limit ourself to r = 0. After summing over k we 
obtain the result 

rn / iS/3h 2iS/3h \ 

g<t>(E) ~ — ttt- ( ~r — 777 + -. 777 J • (58) 

2n \smv/b smt>/3 J 

We next convolve this expression with the Gaussian in order to compare 
with the exact spectrum which has been similarly convolved. To leading 
order in h we obtain for the smoothed "rotational density" 

^ (t?\~ T f sin ( s /3h) „-T 2 w yi8h 2 i sin(2S/3h) 4rV/i87iA fKO \ 

P$,osc\ti) ~ 77— r — ; T-T—e H : — — e . [by 

ZttTi \ smv/b smu/3 J 

The correction to this formula is of order 0(e~ 16T2u,2//18?i2 ) and comes from 
higher repetitions of the orbit. Clearly, with the use of a Gaussian damping 
we have a strong suppression of the longer orbits and higher repetitions of 
short orbits. 

The oscillatory terms may diverge for critical values of the stability angle 
v. The energy dependence of v is shown in Fig. 3. At zero energy it has the 
value of 47r and approaches 3ir near a critical energy close to 0.145 at which 
point the orbit becomes unstable and v becomes complex. The denominator 
of the second term of Eq. ([59]) is zero when v = 3tt so the corresponding con- 
tribution diverges. In this case, we should go to the next order in the saddle 



point approximation as discussed in Ref. |T7|]. However, since the stability 
angle approaches zero abruptly near the critical energy, this divergence is 
only significant for energies very close to the critical energy. 

The final step is to integrate Eq. (|59D with respect to energy to get the 
spectral staircase function. In principle, T, S and v all depend on energy 
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but to a good approximation we can assume S ~ TE and ignore the energy 
dependence of T and v in the integral to obtain the final expression 

ff {g) w _ A I CQS ^) e -^ + cos (25/3^) e _ 4TV/18ft2 \ 
27r \ sinf/6 sin-u/3 / 

In Fig. 4 we show the results of the exact smoothed density of states and 
the results of formulas (53) and (|60"D for two values of the smoothing. For 



w = 0.004, it is clear that the first term of Eq. ( |5"0"D is sufficient to describe the 
oscillations although both terms were included. For w = 0.0015, inclusion of 
the second term is essential. In both cases there is very good agreement over 
the entire energy scale shown. The agreement in Fig. 4 is remarkable given 
the usual inadequacies of the trace formula in mixed systems. We will see 
that for N a (E) the agreement is much poorer. It is worth mentioning that 



the good agreement for N^E) is in contradiction to the result of Ref. [15 
There the authors found that, for unknown reasons, the theory describes 
their numerical results better if the multiplicity of orbit C is taken to be 
three rather than two. Had we done the same, we would have found a fifty 
percent discrepancy in the amplitude of the oscillations rather than the near 
perfect agreement shown. 

The oscillations in the other combination, N a (E), are dominated by the 
unstable orbit B. A similar analysis to that for orbit C gives the result 

N,ME)^-^^e- T2 ^\ (61) 
7rsmh u/4 

where u is the stability of the orbit. The action S and the period T are 
close, but not identical, to those of orbit C. Once again, these are defined 
as the values on the full configuration space. The function u(E) is plotted in 
the bottom of Fig. 3. Eq. (pip only contains the contribution from a single 
repetition. The correction comes from the boundary orbit A and higher 
repetitions of B and is of order 0(e~ T w l 2h ). The result is shown in Fig. 5 
where, unlike before, the agreement is quite poor except in a small region of 
energy between about 0.10 and 0.16. This can be explained by the fact that 
the expression diverges when u is zero. However, as is apparent from Fig. 3, 
u is small for a large range of energy. This is in contrast to the behaviour of 
v which is close to its singular value of 3it only in a narrow range of energy. 
Furthermore, the divergence in Eq. ( |60|) is in the suppressed second term and 
not in the leading term as in Eq. fl6"l"|). 
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7 Conclusion 



In this paper, we have developed a consistent scheme for the evaluation of 
the Weyl density of states for potential problems in the presence of either 
discrete or continuous point group symmetries. For discrete groups, our 
formalism allows us to express the density of states of any irrep as a specific 
combination of contributions from the various classes. By decomposing the 
orthogonal transformations into a product of reflections and rotations, we 
were able to evaluate the contribution coming from any conceivable class of 
group elements. Each class contributes with a relative power of H which is 
simply the codimension of the set of points which is invariant with respect to 
that group element. For an n dimensional system, we are usually interested 
in evaluating at least the first n + 1 terms of the Weyl series so that the 
remainder decreases with h. Because of this, it is important to have a scheme 
to find the h expansion of the contribution from a specific class and also to 
be able to find the contributions from all possible classes. In this paper, we 
have shown how both of these can be done. 

As an example, we showed how these expansions can be used to evaluate 
the Weyl density of states for the two dimensional Henon-Heiles potential 
and excellent agreement was found. In addition, we evaluated the contribu- 
tions from the first periodic orbits and also found good agreement with the 
numerics - again within the symmetry reduced picture. 

The calculation involving continuous rotational groups is subtle because 
of the fact that the limits ft, — > and 6 — > do not commute. For 5*0(2) 
states, we were able to evaluate the contribution of all rotations to the m = 
states and found that these states can be identified exactly with the L = 
manifold in the classical phase space. For other m states, the identification 
is not so clean and we can not evaluate p m {E) exactly but rather must do 
a local average over the index m. We get a similar result for SO (3) and a 
simple extension will give the result for any group of the form SO(n). 

In this work, the discussion was for Hamiltonians of the form H = T + V 
where T depends only on the momenta and V depends only on the coordi- 
nates. This encompasses the majority of physical systems in which one is 
interested. However, one is sometimes interested in more general situations 
where terms involving momenta and coordinates are mixed. One example is 
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a particle in a magnetic field, which has the quantum Hamiltonian 



=i(p- e A(q)) 2 + K(q). (62) 

Another example is the class of systems for which the dynamics is defined 
in terms of group generators, such as the kicked top or various algebraic 



models ||19|| . In these systems, the designation of what is momentum and 
what is coordinate is arbitrary and the Hamiltonian is composed of mixed 
terms. It is possible to extend the results of this paper to such systems. 
However, this is more complicated than what we have done since the Wigner 
transform of the Hamiltonian is no longer found by simply replacing the 
quantum momentum and coordinate operators by their classical counterparts 
as in Eq. (|TT|). Doing so yields only the first term in an ^-expansion of the 
Wigner transform of the Hamiltonian operator. In principle, we would need 
to evaluate all of the terms and keep track of this expansion throughout the 
remainder of the calculation. 

We would like to thank the authors of Ref. [TH] for the use of Fig. 2. In 
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A 2 



1 1 1 
1 1 -1 
2-10 



Table 1: Character table of the group C 3v . The superscripts on the class 
labels indicate the number of group elements in that class. 



Figure 1: The top box shows N^E) and the bottom box shows N a (E) 
as defined in the text. Both refer to the the Henon-Heiles potential. The 
solid lines are the exact numerics while the dashed lines are the analytical 
estimates. Both of these are for h = .004 with a smoothing width w = 0.01. 

Figure 2: The configuration space. The lower left triangle depicts the funda- 
mental domain. The solid lines show the three fundamental periodic orbits 
as they exist in the fundamental domain. The dashed lines show the same 
orbits in the full configuration space. 

Figure 3: The top box shows the period T(E) of the two periodic orbits 
B (dashed line) and C (solid line). The middle section gives the stability 
angle v(E) for periodic orbit C. It varies between 4ir at E = and 3tt at 
E ~ 0.145, as indicated by the dotted lines. The bottom box shows the 
stability exponent of periodic orbit B. 

Figure 4: N^{E) for two smaller values of the smoothing width, w = .004 in 
the upper figure and w = .0015 in the lower figure. The solid lines are the 
exact numerics while the dashed lines are the analytical estimates including 
two repetitions of periodic orbit C. 

Figure 5: N a (E) for a smaller value of the smoothing, w = .0015. The lower 
figure is simply an expansion of the region between E = 0.10 and E = 0.16. 
The solid lines are the exact numerics while the dashed lines are the analytical 
estimates including one repetition of periodic orbit B. 
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